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Abstract. - Tunneling spectra measured on thin epitaxial YBa2Cu307_5 films are found to 
exhibit strong spatial variations, showing U and V-shaped gaps as well as zero bias conductance 
peaks typical of a rf-wave superconductor. A full correspondence is found between the tunneling 
spectra and the surface morphology down to a level of a unit-cell step. Splitting of the zero bias 
conductance peak is seen in optimally-doped and overdoped films, but not in the underdoped 
ones, suggesting that there is no transition to a state of broken time reversal symmetry in the 
underdoped regime. 



The d-wave nature of superconductivity in YBa2Cu307_5 (YBCO) and some other high- Tc 
superconductors is well established by now In particular, various groups have reported 

directional tunneling and point contact spectroscopy measurements that exhibit a dominant 
dx2-y2 symmetry for the order parameter of YBCO ^-||]. However, the spatial variation 
of the electronic properties in the superconducting state and its correlation with the surface 
morphology has not been reported for this system. These properties are reflected in the 
quasi-particle density of states (DOS), which can be measured by tunneling spectroscopy 
[^. Such studies are important for resolving various open problems regarding the nature of 
superconductivity in YBCO. For example, contrary to theoretical predictions zero- 
bias conductance peaks (ZBCP) due to surface Andreev bound state (ABS) were commonly 
observed when tunneling was nominally in the [100] direction [^-||. This was attributed, 
but not yet demonstrated experimentally, to faceting and microscopic roughness that reveal 
the (110) face at which ABS are expected to reside [^||,[l2j. Another intriguing issue is 
that of the emergence of a state with broken time reversal symmetry at the surface due to 
the existence of a subdominant order parameter with 7r/2 phase difference relative to the 



dominant (1^2 _y2 symmetry [|12,13 . This phenomenon manifests itself as a splitting of the 
ZBCP in zero magnetic field at low enough temperatures. Such a splitting was observed by 
some groups but not by others |]6|, ^[l4| and it was shown by Deutscher that the splitting 
might depend on the doping level p, |l5|. Therfore, since the doping level may fluctuate 
spatially over the surface even on a sub- micrometer scale llq, O], the nature of the ZBCP 
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and its splitting should consequently exhibit strong spatial variations. In the present study 
we exploit the high spatial resolution enabled by scanning tunneling microscopy (STM) and 
spectroscopy for the investigation of the above issues. 

Extensive STM studies of high- Tc superconductors have been performed in recent years 
[p^ . These include atomically resolved investigations of the role of disorder in YBCO as 
well as the effects of atomic-size impurities on the local quasi-particle DOS in a single crystal 
of BiSr2CaCu208 (BSSCO) that portray the d-wave nature of this superconductor [p0|,|2l[|. 
However, the effect of larger scale topographic features, such as micro-facets and unit-cell 
steps, on the local DOS was not yet treated, and this is the main focus of our present STM 
study of ultra-thin YBCO films. Our tunneling spectra exhibit pronounced spatial variations 
that are well correlated with the surface morphology. We show that a single unit-cell step 
on the a&-surface is sufficient to induce a ZBCP when tunneling along the c axis. Moreover, 
splitting of the ZBCP was found in various locations (but not everywhere) on nominally (as 
prepared) overdoped and optimally-doped samples, but was never observed in the underdoped 
films, consistent with ref. |]l5|] . 

Epitaxial YBCO films of 200 - 300 A thickness were grown on (100) SrTiOs wafers by 
laser ablation deposition, with c-axis orientation normal to the substrate. Films of different 
doping levels were obtained by the use of different annealing conditions immediately after the 
deposition process. Annealing at 0.8 atm oxygen pressure with a dwell of Ih at 450 C yielded 
optimally or slightly overdoped films with ^ 90 K, while a cool-down under 20 niTorr 
oxygen flow at 900 C/h yielded an underdoped phase of ~ 60 — 70 K. The of our films 
was uniform up to about 0.5 K, as measured at 10 different points over the 1 crri^ area of 
the wafer. X-ray diffraction analysis and the structure of the R(T) curves (not shown here) 
confirmed the doping levels of these phases. At a relatively low deposition rate and small laser 
power, the films consisted of crystallites having well defined facets, as depicted in fig. [l^. The 
crystallite structure became less clear upon increasing of the laser power, and the dominant 
surface structures revealed by the STM images were voids and island of one unit-cell height 
(fig. |a). The rich surface morphology of these films enabled the observation of a wealth of 
spectral structures and the study of their spatial dependence, as shown below. The films were 
taken out from the growth chamber and transferred in a dry atmosphere to a cryogenic STM 
within a few hours. The onset of the superconducting transition did not change during this 
transfer process, but the transition width increased considerably from about 1 K for freshly 
prepared films to around 10 K after completing the STM measurements. The tunneling 
spectra (dl/dV vs V characteristics) were acquired either directly by the use of conventional 
lock-in technique, or by numerical differentiation of the measured I-V curves, with similar 
results obtained by both methods. We have checked the dependence of the tunneling spectra 
on the voltage and current setting (i.e., the tip-sample distance) and found it not to affect 
the measured gap and ZBCP features. However, the setting did influence the background, 
in particular the degree of asymmetry the spectra exhibit, as noted previously by various 
groups [?]. 

A typical correlation between the tunneling spectra and the local surface morphology is 
depicted in fig. |l|, measured on an optimally-doped film, which focuses on a nearly square 
shaped crystallite about 50 nm long and 17 nm high (fig. |l]a). On top of this crystallite, 
V-shaped gap structures were observed, as shown by the curve acquired at position labeled 
1 in this image, presented in fig. |^b (solid line). In general, spectra measured on such 
positions showed fiuctuations in the degree of asymmetry, in the sharpness of the gap-edge 
shoulders and in the possible existence of small zero-bias anomalies (see below), but not in 
the general V-shaped structure as expected for tunneling in the c direction. This is due to 
the finite tunneling-cone width in the STM experiment that allows coupling to the electronic 





Fig. 1 - Correlation between tunneling spectra (at 4.2 K) and the surface morphology of a thin 
YBCO film, (a) 50a;50 nm? STM topographic image showing four crystallites and the interfaces 
between them. The numbers 1-3 denote the positions where spectra b, d and f were acquired, (b- 
f) Tunneling spectra measured at different locations (black lines), where (b-e) were taken on the 
same sample and (f) on a different one. The dotted curves are fits to the theory of tunneling into 
a d-wave superconductor with order parameters, A^j, denoted in the figures. The spectra exhibit 
various structures that manifest rf-wave superconductivity: (b) V-shaped gap for the curve measured 
on top of a crystallite (position 1), typical for c-axis tunneling, (c) Asymmetric V-shaped gap with 
a small zero-bias structure (measured on top of another crystallite) ; the 1- V curve is presented in the 
inset, (d) U-shaped gap for the spectrum taken at the interface position 2, where the (100) face is 
revealed, (e) ZBCP inside a gap structure for the curve measured near a corner position 3, showing 
contributions of (100), (110) and (001) faces, (f) A pure ZBCP structure taken on a (110) facet. 



states in the Cu02 planes. The general structure of this dl/dV-V curve is well reproduced 
by a simulation based on the formalism for tunneling into a c?-wave superconductor [pO|,pT[|, 
assuming tunneling in the c direction (dotted curve). In many cases a small structure (far 
from being a well-developed ZBCP) appears at around zero bias, as exhibited by the curve in 
fig. I^c, taken on top of a different crystallite. Similar structures were found by other groups 
on c-axis BSSCO crystals and attributed to impurity-induced quasi-particle excitation at the 
Fermi level [|o[|l|,|2j,^. Enhanced gap-edge shoulders, as seen in the spectrum, have also 
been observed previously j|,|l^,|2^ and were related to a van-Hove singularity in the electronic 
DOS 

The tunneling spectra obtained when positioning the tip above a crystallite edge look 
considerably different when either the (100) or the (110) planes are probed. A nearly U- 
shaped gap structure is obtained at position 2 (fig. ^), as well as in any other position 
along the same facet, suggesting that the (100) face is exposed at this edge, as confirmed by 
the relatively good fit that takes into account such a contribution. The spectra obtained in 
position 3 show a clear ZBCP inside a gap (fig. |^). The ZBCP is attributed to ABS formed 
at the exposed (110) facet. Typically, however, the tunneling spectra acquired on larger (110) 
crystalline edges exhibited more pronounced ZBCP structures, as depicted in fig. The 
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Fig. 2 - The effect of steps of unit-cell height on the tunneling spectra measured on a slightly overdoped 
film, (a and b) Topographic images of size 120a;220 and i20a;i20 nm^ , respectively, showing monocell 
steps, islands and voids on the (001) surface, (c and d) Tunneling spectra (at 4.2 K) taken on 
different locations as denoted in image b. Curve 1, measured far from the step, has a V-shaped gap, 
manifesting c-axis tunneling. Curves 2 and 3, measured near the steps, exhibit ZBCPs of different 
magnitudes inside U and V-shaped gaps, respectively, (e) Two highly asymmetric tunneling spectra 
showing pronounced ZBCPs, taken at a different step on the same sample. To demonstrate the spatial 
variation of the spectra within this scan range we plot in the insets of (c) and (d) additional data 
acquired at analogous positions as those presented in the main frames. 



fit to the curve in fig. |i|e was obtained by taking contributions of (100), (110) and (001) 
faces, consistent with the corner position (3 in frame a), while that in fig. has only a (110) 
component. 

The main topographic features observed on films that were deposited using a higher laser 
power were islands and voids of unit-cell height on the (001) surface, as shown in image 
Image ||b focuses on a corner between two steps, each of a single unit-cell height, around 
which the spectra presented in figures |2[;-|^ were measured. Generally, tunneling spectra 
taken far enough (a few nm) from a step or a surface defect (such as holes of sub unit-cell 
depth appearing in the image as black dots) exhibited a V-shaped structure, such as curve 
1 in fig. In the vicinity of steps and "point-defects", ZBCPs of different magnitudes 

embedded in relatively symmetric gap structures were observed. This is exhibited by curve 2 
(fig. Il^), taken near a (100) mono-step facet and showing a U-shaped gap, and by curve 3 (fig. 
gd), acquired on a (110) mono-step facet, portraying a ZBCP within a V-shaped gap. The 
difference in the apparent gap between curves 1 and 2 may be due either to local variations 
in the oxygen content or to a life-time broadening effect. Such variations were observed also 
in STM measurements on single crystals Q. To demonstrate the typical variance of the 
spectra within the area shown in image |^(b) we plot in the insets of (c) and (d) additional 
characteristics acquired at analogous positions to those presented in the main frames. It is 
evident that the overall shapes are reproducible along a specific topographic feature. To the 
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Fig. 3 - Splitting of the ZBCP. (a) Two tunneling spectra measured at the same position, one at 
4.2 K (solid line) and the other at 10 K (dotted). The clear split seen at 4.2 K is completely washed 
out at 10 K, signifying a transition into a broken time reversal symmetry state, (b) Three tunneling 
spectra showing (1) the largest splitting we have detected (3.2 meV peak to peak), (2) the smallest 
splitting (1.6 meV) and (3) a pronounced asymmetric curve. The dotted line is a fit assuming a state 
of broken time reversal symmetry having a d^2_y2 + is order parameter, with Ad — 30 meV and 
A, = 0.04Ad. 



best of our knowledge, this is the first time that a clear correlation between tunneling spectra 
and unit-cell step structure and orientation is demonstrated. This indicates that even a single 
step of one unit-cell height can alter considerably the measured local quasi-particle excitation, 
pointing out to the apparent discrepancy between theory and experiment, discussed above. 

We have also observed tunneling spectra that deviate from the standard d-wave behavior. 
For example, highly asymmetric ZBCP structures were found on the same sample in similarly 
looking areas (fig. |^) . Such spectra were predicted to result from scattering by atomic-scale 
impurities |^,|2^. Unfortunately, we can not yet specify any typical defect structure, or any 
special topographic feature, which can be directly related to the appearance of these highly 
asymmetric tunneling spectra. 

We now turn to the issue of splitting of the ZBCP, a manifestation of a state with broken 
time reversal symmetry, which we have observed only in overdoped and optimally doped films. 
In fig. ^ we present two tunneling spectra obtained at the same position on a nominally 
overdoped film, one taken at 4.2 K, exhibiting a clear splitting of the ZBCP (~ 2.6 meV 
peak to peak), and the other at 10 K, where the split is completely washed out. (The 
lateral position of spectra acquisition was determined at each temperature via the topographic 
images, although slightly different tip-sample separations may have yielded a small difference 
in the background.) This temperature dependence indicates that the splitting is due to a 
broken time reversal symmetry phase. The peak to peak separation we have found in all 
of our measurements ranged between 1.6 meV (curve 2 in fig. ^a) to 3.2 meV (curve 1). 
The theoretical fit to curve 1 was obtained using a d^2_y2 + is order parameter [p^ , where 
the magnitude of the subdominant s component, Ag, is 4% of the dominant one, A^. We 
note that we cannot clearly distinguish between the is ||l^ and the idxy (l^ scenarios for 
the subdominant order parameter, although recent experiments seem to support the latter 
case 

Interestingly, the splitting into two well resolved peaks we typically observed at 4.2 K was, 
in many cases, more pronounced than the splitting into two weakly resolved peaks measured in 
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macroscopic planer junction at zero field The reason for this is that in the macroscopic 

tunneling junctions the dl/dV-V curves display an average over large areas, whereas the STM 
allows probing of specific locations where the splitting is more pronounced. Local probing 
enables also the observation of highly asymmetric spectra that exhibit a state of broken time 
reversal symmetry, such as curve 3. 

The fact that splitting was not observed in underdoped samples indicates that sponta- 
neous time reversal symmetry breaking may take place only beyond some critical doping 
level, as pointed out by Deutscher jisj. This may be the reason for the contradictory re- 
ports of various groups regarding the observation of splitting in zero magnetic field, discussed 
above. The spatial variations in magnitude of the splitting seen in our nominally overdoped 
and optimally doped samples may thus be due to fluctuations in the local hole doping. The 
wide superconducting transition exhibited by our samples after completing the STM measure- 
ments, as compared to the freshly prepared films, indeed indicates that the samples became 
non-homogeneous electronically, probably due to inhomogeneous oxygen distribution. Clear 
evidence for spatial fluctuation in the hole doping was previously observed by STM in other 
high-Tc superconducting systems [|l6|,|l7||. 

In summary, tunneling spectroscopy measurements on ultra-thin epitaxial YBCO films 
reveal a rich variety of rf-wave spectral structures that show pronounced spatial variations 
with nearly one-to-one correspondence to the surface morphology. Even a single unit-cell step 
is found to significantly affect the local spectra. The doping and position dependence of the 
ZBCP splitting suggest that doping plays a major role in the transition to a state with broken 
time reversal symmetry, which can occur only beyond some doping level close to optimal 
doping. 

* * * 
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